Astronomy & Astrophysics manuscript no. msmith'submitfinal 


© ESO 2012 


February 1, 2012 





The relationship between 7 Cassiopeiae's X-ray emission and its 

circumstellar environment 

M. A. Smith 1 , R. Lopes de Oliveira 2 ' 3 , C. Motch 4 , G. W. Henry 5 , N. D. Richardson 6 , K. S. Bjorkman 7 , Ph. 
Stee 8 , D. Mourard 8 , J. D. Monnier 9 , X. Che 9 , R. Biicke 10 , E. Pollmann 11 , D. R. Gies 6 , G. H., Schaefer 6 , T. 
ten Brummelaar 6 , H. A. McAlister 6 , N. H. Turner 6 , J. Sturmann 6 , L. Sturmann 6 , and S. T. Ridgway 12 

1 Catholic University of America, 3700 San Martin Dr., Baltimore, MD 21218 USA; e-mail: msmith@stsci.edu, 

2 Universidade Federal de Sergipe, Departamento de Fi'sica, Av. Marechal Rondon s/n, 49100-000 Sao Cristovao, SE, 
Brazil 

3 Universidade de Sao Paulo, Instituto de Fi'sica de Sao Carlos, Caixa Postal 369, 13560-970, Sao Carlos, SP, Brazil 

4 Universite de Strasbourg, CNRS, Observatoire Astronomique, 11 rue de l'Universite, 6000 Strasbourg, France 

5 Center of Excellence in Information Systems, Tennessee State University, 3500 John Merritt Blvd., Nashville, TN 

6 Center for High Angular Resolution Astronomy, Department of Physics and Astronomy, Georgia State University, 
P.O. Box 4106, Atlanta, GA 30202-4106, USA 

7 Ritter Astrophysical Research Center, Department of Physics & Astronomy, University of Toledo, 2801 W. Bancroft, 
Toledo, OH 43606, USA 

8 Laboratoire Lagrange, UMR 7293 UNS-CNRS-OCA, Blvd. l'Observatorie, B. P. 4229F, 06304 Nice Cedex 4, France 

9 Department of Astronomy, University of Michigan, Ann Arbor, MI 48109, USA 
10 Anna-von-Gierke-Ring 147, 21035 Hamburg, Germany 

Emil-Nolde-Str. 12, 51375 Leverkusen, Germany 
12 National Optical Astronomical Observatory, 950 North Cherry Ave., Tucson, AZ 85719, USA 

Received October dd, 2011; accepted ????? DD, YYYY 

ABSTRACT 

7 Cas is the prototypical classical Be star and is recently best known for its variable hard X-ray emission. To elucidate 
the reasons for this emission, we mounted a multiwavelength campaign in 2010 centered around four XMM-Newton 
observations. The observational techniques included long baseline optical interferometry (LBOI) from two instruments 
at CHARA, photometry carried out by an Automated Photometric Telescope and Ha observations. Because 7 Cas is 
also known to be in a binary, we measured radial velocities from the Ha line and redetermined its period as 203.55±0.20 
days and its eccentricity as near zero. The LBOI observations suggest that the star's decretion disk was axisymmetric 
in 2010, has an system inclination angle near 45°, and a larger radius than previously reported. In addition, the Be 
star began an "outburst" at the beginning of our campaign, made visible by a brightening and reddening of the disk 
during our campaign and beyond. Our analyses of the new high resolution spectra disclosed many attributes also found 
from spectra obtained in 2001 (Chandra) and 2004 (XMM-Newton). As well as a dominant hot (~ 14keV) thermal 
component, the familiar ones included: (i) a fluorescent feature of Fe K even stronger than observed at previous times, 
(ii) strong lines of NVII and NeXI lines indicative of overabundances, and (iii) a subsolar Fe abundance from K- 
shell lines but a solar abundance from L-shell ions. We also found that two absorption columns are required to fit the 
continuum. While the first one maintained its historical average of lx 10 21 cm -2 , the second was very large and doubled 
to 7.4xl0 23 cm -2 during our X-ray observations. Although we found no clear relation between this column density 
and orbital phase, it correlates well with the disk brightening and reddening both in the 2010 and earlier observations. 
Thus, the inference from this study is that much (perhaps all?) of the X-ray emission from this source originates behind 
matter ejected by 7 Cas into our line of sight. 

Key words, stars: individual: 7 Cas stars: emission-line, Be; X-rays: stars, stars: activity, stars: circumstellar matter, 



1. Introduction 

Discovered in 1866 by Secchi from the emission in its Ha 
spectral line, 7 Cas became the first member of the now well 
known group of "classical Be" variables (Porter & Rivinius 
2003). These stars are also known for their rapid rotation, 
continuum excess in the infrared, and by the presence of 
strong, blueshifted, and often variable, absorption compo- 
nents in their UV resonance lines. The latter were stud- 
ied extensively in 7 Cas and other Be stars at times inten- 
sively during the 1978-1996 operation of the International 
Ultraviolet Explorer (IUE). Fortunately, soon after of the 



termination of the operations of this satellite, other tech- 
niques were being introduced to astronomy and have ad- 
vanced the study of the Be stars and their ejected Keplerian 
disks. These developments include the operation of Long 
Baseline Optical Interferometry (LBOI) facilities, the de- 
velopment of automated photometric (robotic) telescopes, 
and an increasingly relevant contribution in spectroscopy 
by amateur telescopes with their own local facilities. The 
present study makes use of these elements. 

Through the last decade astronomers have mounted a 
series of multi-wavelength campaigns on the prototype of 
the "classical" Be stars, 7 Cas (HD5394, type B0.5IVe). 
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Even more than most classical Be stars, 7Cas is vari- 
able over a number of wavelengths and timescales - see 
Harmanec (2002). The study of 7 Cas developed renewed fo- 
cus by the discovery that it emits X-rays with unusual prop- 
erties and in particular by White et al.'s (1982) suggestion 
that these emissions arise from infall onto a degenerate com- 
panion. We now know that this star is in fact the primary in 
a nearly circular, 203-day binary system (Harmanec et al. 
2000, Miroshnichenko et al. 2002). It may be important that 
its spectral lines are highly rotationally broadened (vsin i 
400-441 kins" 1 ; Harmanec 2002, Frerhat et al. 2005). With 
this broadening and a probable rotational period of 1.2158 
days (Smith, Henry, & Vishniac 2006; "SHV"), this star is 
rotating at nearly its critical rate. 

7 Cas is also the prototype of a group that so far includes 
eight other Galactic "7 Cas X-ray analogs." These include: 
HD 110432 (Smith & Balona 2006; "SB06"), SS397, "Star 
9" in NGC 6649 (USNO 0750-13549725; Motch et al. 2007), 
HD 161103 & SA049725 (Lopes de Oliveira et al. 2006; 
"L06"), HD 119682 (Rakowski et al. 2006, Safi-Harb et al. 
2006), XGPS-36 (Motch et al. 2010), and HD 157832 (Lopes 
de Oliveira & Motch 2011). The X-ray properties of this 
group are not shared by the other classical Be stars. 

In the years since the White et al. study, much effort 
has been devoted to characterizing the X-ray attributes of 
7 Cas itself because of its brightness. We know now that 
these properties include: (1) a rather high mean L^/L ~ 
10~ 6 (but at least 10 times lower than Be X-ray binaries), 
(2) a light curve variable on a number of timescales, and (3) 
a thermal spectrum consisting of 3-4 distinct plasma com- 
ponents. It is dominated by a hot component with temper- 
ature kThot of 12-14 keV and subject to absorption by two 
independent columns of cold (ambient) gas (Smith et al. 
2004, Lopes de Oliveira 2010). It should also be remarked 
that the X-ray emission consists of an underlying "basal" 
component, which varies on a timescale of hours, punctu- 
ated by rapid flares that according to plasma cooling argu- 
ments alone must be produced in a high density environ- 
ment (Smith, Robinson, & Corbet 1998, "SRC")EI 

The X-ray light curve also exhibits variations on at 
least three timescales. Much of the behavior of the X-ray 
light curve of 7 Cas is now well investigated by a series 
of relatively homogeneous RXTE studies (SRC, Smith & 
Robinson 1999, Robinson & Smith 2000, Robinson, Smith, 
& Henry (2006; "RSH"), and SHV). A rapid variability, 
first described in detail by Murakami et al. (1986), takes 
the form of ubiquitous flares on timescales of 2-3 minutes 
down to 4 seconds (the instrumental detection threshold). 
These flares are "mild," only seldom having amplitudes of 
more than 2 x above the background "basal" flux. The basal 
fluxes undergo undulations on timescales as short as 20 min- 
utes but more typically a few hours. Thirdly, a long X-ray 
cycle of ^70 day observed by RSH and SHV seem to be 
well correlated with long periods at optical wavelengths. 
The total range of these cycles in the optical range is about 
50-91 days. 

In this paper we seek to understand the X-ray emission 
generation in terms of its possible relationship to the gen- 
eral Be circumstellar environment, including both the de- 
cretion disk and its widely spaced binary companion. The 



1 We define a flare as a sudden and short-lived increase in 
X-ray flux. This term does not necessarily connote a magnetic 
origin for them as on the Sun. 



impetus for this study was the different amounts of soft 
X-ray flux in earlier high-resolution spectra by Chandra 
and the XMM-Newton satellites (Smith et al. 2004, Lopes 
de Oliveira et al. 2010). The different soft X-ray flux lev- 
els are due to a density range (by a factor of 40) in one 
of two absorption columns used in our models to fit the 
full XMM-Newton spectrum (the other column length solu- 
tion was stable). We noticed that the 2001 Chandra obser- 
vations were obtained during the binary phase of inferior 
conjunction, that is with the Be star passing through the 
foreground in its orbit. To see if the column density de- 
pends on binary phase, we requested four new exposures in 
2010 for inferior conjunction phase and extending almost 
to quadrature. The question we posed was whether extra 
absorption would recur at this phase, e.g. due to absorption 
by matter near the Be star. 

2. Description of the circumstellar matter 

Many of the optical and X-ray variations from 7 Cas are 
likely to be caused somehow by interactions of matter 
released by the star with its disk or with previous ejecta. 
Since the present study sheds new light on some these 
variations we should first review the status of the circum- 
stellar (CS) structures. 

Geometry of the disk: 

For several decades the emission from the decretion disk 
of 7 Cas has varied irregularly, presumably because of dis- 
continuous feeding of matter to its decretion disk (Doazan 
1982, Doazan et al. 1987). The emission in the continuum 
is formed by free-free and bound-free process and in the hy- 
drogen lines by recombination. Stee & Bittar (2001) found 
that the disk is optically thin at visible wavelengths but op- 
tically thick in the near infrared and in the Ha line out to 
greater disk radii. Their results indicated that the disk con- 
tributes ;> 10% to the total flux of the star-disk system in 
the Johnson V band but only 3-4% in B. Spectroscopic and 
interferometric arguments suggest that the disk and rota- 
tional axes coincide, and both techniques indicate that this 
angle, the inclination, has an intermediate value. For ex- 
ample, in their initial interferometric work in Ha, Mourard 
et al. (1989) determined a disk inclination angle of ~45° 
from its measured ratio of the minor to major axis on the 
sky. Observations at this wavelength by Quirrenbach et al. 
(1997) and Tycner et al. 2006), with the Mklll and NPOI 
interferometers, respectively, as well as in the 2.1 /xm K' 
continuum band by Gies et al. (2007) have determined in- 
clinations in a range i = 46-55°. 

LBOI studies have also determined the disk size of 
7 Cas. Tycner et al. Ha found a diameter of the ma- 
jor axis of 3.67 mas. This compared well with the 
Quirrenbach et al. (1997) measurement of 3.47 mas. 
In this paper we will adopt a uniform angular di- 
ameter for the Be star of 0.44 mas. Therefore, 
these measurements correspond to about 8R». In 
addition, ost kinematic studies of Be disks have found that 
they have Keplerian rotation (e.g., Meilland et al. 2007). 
Mourard et al. (1989) could show that the 7 Cas disk 
rotates but were unable to specify its velocity-distance 
relation. 

Nonsecular variations imposed by the disk: 

Historically 7 Cas has been among the \ of classical Be 
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stars for which the violet ("V") and red ("i?") compo- 
nents of their Balmer emission lines exhibit cyclic variations 
with timescales of years (e.g., Hummel 2000). These cycles 
are believed to be caused by the precession of a 1-armed 
density structure from the excitation of nonaxisymmctric 
disk modes (Okazaki 1991). They were first observed in 
Balmer lines of 7 Cas in 1969 (Cowley, Rogers, & Hutchings 
1976). Continued observations in the 1980s-1990s disclosed 
a strong cycle with a length of 5-6 years. Berio et al. (1999) 
traced the precession of this arm over a few years by LBOI 
observations in Ha. Rittcr Observatory spectra obtained by 
KSB show that the cycle was still present in 2000. Later ob- 
servations by Rivinius (2007) indicate that the cycle weak- 
ened through 2003. By 2005 the V/R variations recorded by 
Ondrejov Observatory spectra were "very weak if present 
at all" (Nemravova 2011). Spectra obtained through the 
Be Star database ("BeSS", Neiner et al. 2011) from 2009 
to the present, including the epoch of our X-ray observa- 
tions, show no variations either. 

Typical of other classical Be stars observed at low or 
intermediate inclinations, the UV resonance "wind" lines 
in the 7 Cas spectrum indicate the presence of a variable 
high velocity wind. Henrichs, Hammerschlag-Hensberge, & 
Howarth (1983) found that optical depths in the wind lines 
can vary by a factor of up to 20 over time. The optical 
depths in the wind components increase by a mean factor 
of 2-3 times from V/R <1 to V/R >1 phases. Telting & 
Kaper (1994) advanced a picture in which the wind ablates 
disk matter, thereby focusing the wind and increasing col- 
umn absorption at "blue" {V/R > 1) phases. Beyond this 
modulation, changes in wind line strength can vary by up 
to a factor of five in a few weeks. 

The cause of changes in the wind over these shorter 
timescales is a still unsolved mystery and is probably re- 
lated to how Be stars change their local environment. In a 
distinctly different scenario from the Telting-Kaper picture, 
Brown et al. (2008) have suggested that a global magnetic 
field from the star focuses the wind towards the disk, thus 
producing a compressed and "magnetically torqued disk." 
In such a picture one can imagine that time-dependent 
changes in the field geometry would focus the wind to vary- 
ing degrees across our line of sight at different times, pro- 
ducing variations in the UV wind absorption components. 
This mechanism could not work for 7 Cas globally because 
a dipolar field would produce unobserved strong periodic 
emission and absorption variations over the rotation cycle, 
as observed only in Bp stars. Still, it is conceivable that 
magnetic loops from chaotic surface fields could divert wind 
flows on smaller scales. 

Ha emission and optical brightening as well as red- 
dening episodes in the optical records of 7 Cas and other 
Be stars occur over timescales as short as a few weeks. In 
addition, Carciofi et al. (2012) have simulated a ^100 day 
dissipation of the disk of the Be star 28 CMa with a viscous 
disk model and concluded the high Shakura-Sunyaev 
viscosity implied, a = 1.0±0.2, requires that the mass 
injection rate into this disk exceeds the wind mass loss 
rate by about an order of magnitude. These considerations 
suggest that an unknown deus ex machina operates 
at irregular intervals that causes discrete ejection events 
("outbursts") in at least some Be stars. In ij7]we will return 
to this topic. For now we will loosely use the term "wind" 
to indicate any type of ejection mechanism, although they 



may well involve more than radiative and centrifugal forces. 
Corotation and migrating subfeatures: 

A number of observers have noted the existence of nar- 
row absorptions called "migrating subfeatures" (msf) that 
run blue to red across optical and UV lines of the 7 Cas 
spectrum. (Yang et al. 1988, Horaguchi et al. 1994, Smith 
1995, Smith & Robinson 1999). In optical lines msf do not 
generally appear during every night of observation. When 
they do appear a new feature moves across a given point 
on the profile every 4-5 hours and at an acceleration rate 
of +95 kms _1 hr _1 . The msf have also been observed in 
spectral lines of the 7 Cas analog HD 110432 (SB06) and 
the magnetically active, rapidly rotating K0 star AB Dor 
(Collier Cameron & Robinson 1989) and are ascribed to 
forced-corotation of dense clouds over the star's surface by 
magnetic fields. Indeed, Donati et al. (1999) have confirmed 
the detection of a complex field topology on the surface of 
AB Dor, reinforcing the corotation interpretation. 

The character of the spectroscopic msf makes it un- 
likely that they arise from nonradial pulsations (NRP). 
Recent photometric monitorings of Be stars have disclosed 
excited modes for which the periods tend to be clumped 
into noncoherent groups. For early-type Be stars they gen- 
erally are excited with longer periods (~0.5 day -1 ) than 
the times between passages of successive jns/absorptionsQ 
The observed NRP modes are attributed mainly to rota- 
tional splittings of g-modes with |m| « 2 (e.g., Walker et al. 
2005, Dziembowski et al. 2007, Balona et al. 2011) or to in- 
homogeneities of orbiting disk matter (Balona 2009). Both 
their short lifetimes and their appearance as high frequency 
"wiggles" in the line profiles make the msf inconsistent with 
low-|m| modes. 

The interpretation of forced corotation of clouds was 
also advanced to explain dips in UV continuum light curves 
of 7 Cas (Smith, Robinson, & Hatzes 1998; "SRH"). The 
span of these dips is too brief to be caused by advection 
of spots across the star's surface. They can be explained, 
however, by the presence of translucent clouds corotating 
within a few tenths of the stellar radius of the surface. Also, 
the appearance of UV light curve dips correlates well with 
the increases in X-ray flux over several hours, according to 
a simultaneous X-ray /UV campaign described by SRC. 

3. Origin of the X-ray emissions 

The ultimate goal of our ongoing studies of this topic is to 
define the mechanism responsible for the unique emission 
of the hard and secondary soft X-ray emissions that set the 
7 Cas variables apart from the emissions of other high-mass 
stars. Most recent suggestions for this origin have centered 
around two competing ideas. 

According to the first idea the hard X-rays result from 
the conversion of gravitational energy from accretion of the 
Be star's disk onto the surface of a degenerate companion 
(white dwarf or neutron star). The primary rationale for 
this idea is the resounding success in explaining X-ray ac- 
tive Be stars in the Galaxy and Magellanic Clouds. Indeed, 

2 Exceptions to this statement are found in the spectral vari- 
ations of n Cen and 8 Sco, which can exhibit modes with pe- 
riods in this range in addition to others (Rivinius et al. 2003, 
Smith 1986). In these cases the appearance of broad sinuous 
NRP bumps on line profiles is very different than narrow msf 
depicted in the references we have cited. 
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most of these X-ray sources are known to be binaries be- 
cause of the detection of X-ray pulses due to a neutron star. 
In the case of the 7Cas stars no such pulses are known, 
and the direct detection of binarity through radial velocity 
variations has so far been confined to 7 Cas itself. As al- 
ready indicated, the near-circular orbit of this system hints 
that the secondary is not a neutron star. In this connec- 
tion it may be significant that two members of the 7 Cas 
group, "Star 9" of NGC6649 and HD 119682 seem to be 
blue stragglers of clusters having ages of about 50 Myr and 
40Myr, respectively (Marco, Negueruela, & Motch 2007, 
Marco et al. 2009). Moreover, HD 110432 may be a mem- 
ber of NGC 4609 (age of 60 Myr; Feinstein & Marraco 1979) , 
and if so would be a straggler. Thus, while unproven, it is 
conceivable that all 7 Cas stars are blue stragglers. This cir- 
cumstance could mean either that the current degenerate 
secondary star is the site of the X-ray emission or that it 
simply provides a mechanism to spin up the outer layers of 
the Be star to near criticality, thereby somehow triggering 
a mechanism for the production of the X-ray emissions. 

Another argument for tying the X-ray emission to white 
dwarfs is the multi-thermal nature of the X-ray plasma ob- 
served in some cataclysmic variables, particularly the mag- 
netic "polars" (e.g., Mukai et al. 2003). However, upon close 
inspection there are also significant differences between the 
X-ray characteristics of these objects and 7 Cas. These in- 
clude in the former case a general lack of ubiquitous flaring 
confined to the hot X-ray regions and their exhibiting no 
differences in the Fe abundance determined from L and 
K shell ions (e.g., Itoh et al. 2006). In addition, even in 
their high states the hard components of polars have lower 
fluxes, L x -10 30 -10 32 ergs' 1 (e.g., Watson k King 1987) 
than 7 Cas . If accreting white dwarfs are the source of the 
X-rays, it is likely that they would represent a new subclass 
of such objects. 

A very different concept to explain the X-ray flux of 
7 Cas is called the "magnetic star-disk interaction" model 
and was proposed by RSH. These authors pointed to the X- 
ray flares noted above, the correlation of rapid variations of 
X-ray flux and UV and optical diagnostics associated with 
matter in the proximity of a Be star, some of which were 
noted above, as well as the correlation of the long X-ray 
and optical cycles. In this picture magnetic lines of force 
from active regions on the rapidly rotating Be star and its 
slower rotating disk occasionally get entangled. The result 
is that they stretch, sever, and upon reconnection fling high 
energy particles in various directions including toward the 
Be star. The impacts of these particles on random points 
on its surface heat photospheric plasma, and they become 
visible as "flares" - though not flares in the solar sense of re- 
leasing local magnetic energy. This matter expands to lower 
density regions above the flare site and continues to glow 
as "basal" emission. The observation of X-ray flux "lulls" 
over several epochs suggests that this does not happen per- 
fectly continuously but instead can be subject to relaxation 
cycles. 

The heart of the evidence for the Be star-disk magnetic 
interaction picture lies in observations of rapid correlations 
between X-ray flux and fluxes in other wavelength regimes 
and, secondarily, on inferences of magnetic activity on the 
star's surface. The first well documented example of such a 
correlation was the report by Slettebak & Snow (1978) of 
an Ha emission "flare" that occurred at the same time as 
short-lived occurrences of emissions in scanned components 



of the Mg II and Si IV lines by the Copernicus UV-scanning 
spectroscopic satellite. This could not be dismissed as a 
(double) instrumental "fluke" because Copernicus, using its 
coaligned X-ray telescope, also observed an X-ray during 
this time (Peters 1982). On another occasion, Smith (1995) 
observed a brief "emission flare" in the He I 6678 A line, 
although no other wavelength regions were monitored. 

To pursue this evidence SRC mounted a 21 hour 
simultaneous campaign using the Rossi X-ray Timing 
Explorer (RXTE) and the Hubble Goddard High Resolution 
Spectrograph in 1996 March. They identified several correla- 
tions between X-ray and UV continuum flux as well as var- 
ious spectral lines expected to be formed in the neighbor- 
hood of a B0. 5 star (e.g., Smith & Robinson 2003). These 
variations are the strongest single set of arguments that the 
X-rays are excited near the surface of the Be star. 

The correct explanation for the production of the hard 
X-rays will have interesting ramifications for our under- 
standing evolutionary and possibly the angular momentum 
history of the X-ray emitting star. The elucidation of this 
true mechanism will require the pulling together of diverse 
properties of the Be star, its binary orbit, and circumstel- 
lar matter from observations at several wavelength regimes 
over as contemporaneous a time interval as possible. The 
following describes a comprehensive effort to this end. 

4. Observations 

We first outline the observational techniques used in this 
paper. We begin by summarizing the LBOI results and dis- 
cussing results from automated photometry and Ha line 
spectroscopy. 

4.1. The CHARA Long Baseline Optical Interferometer 

Our LBOI observations were conducted in the infrared us- 
ing Georgia State University's Center for High Angular 
Resolution Astronomy (CHARA) array facility on Mt. 
Wilson. CHARA is an interferometric array of six 1-meter 
telescopes configured in a Y-shaped pattern with base- 
lines of up to 330 m for maximum spatial resolution. Our 
observations utilized two new-generation instruments for 
imaging of 7Cas's decretion disk: the Michigan Infrared 
Combiner (MIRC) in the H-band (1.64 /im) and the Visible 
spEctroGraph and polArimeter (VEGA), which became 
operational in 2009 (Mourard et al. 2009). VEGA obser- 
vations were conducted in Ha and also near this line in 
~200A regions centered at 6460 A, 6690 A, and 6728 A. In 
the first case the results from the three wavelength settings 
were averaged together and we call the result the "R"-band 
measures. The Ha line fluxes were determined by stepping 
in 3 A steps across the 6550-6575 A region. The observa- 
tions were executed in order to describe the geometry of the 
Be disk. The VEGA kinematical program was conducted to 
test the disk rotation-radius relation. These observations 
are summarized in the following and are detailed in Stee et 
al. (2012; "Paper 2"). 

The MIRC instrument was used in a 4-telescope ar- 
rangement, resulting in the simultaneous measurement of 6 
visibilities, 4 closure phases and 4 triple amplitudes. These 
observations used a broadband H filter along with a prism 
to disperse the light into 8 spectral channels (R~45). This 
light is refocused to produce interference fringes received by 
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a "PICNIC" camera designed to maximize rapid data read- 
out and minimize fringe decoherence due to atmospheric 
seeing. Further details on this electronic and readout sys- 
tem may be found in Monnier et al. (2004, 2008), Pedretti 
et al. (2009), and Che et al. (2010). The maximum angular 
resolution in the H band using the MIRC at CHARA 
is 0.4 milliarcsec, or coincidentally about the expected an- 
gular size of 7 Cas, the Be star. The observations for this 
study were carried out on 2010 Aug. 5, Aug. 8, Sept. 24, 
and Nov. 4-5. This range of dates includes the range of 
the XMM-Newton OBS 3 and 4 referred to below and the 
next two months. These observations were followed almost 
immediately by VEGA observations. The MIRC baselines 
ranged from 34 m (S1S2) to the maximum 330 m (S1E1). 

The VEGA instrument is different from MIRC in work- 
ing with visible wavelength (4500-8500 A) light at "low, 
medium, or high" dispersions (1,700, 5,000, 30,000) but en- 
ables the combining of light from up to four CHARA tele- 
scopes and measurements of visibility modulus, phases, and 
phase closures as a function of baseline and wavelength. Its 
angular resolution is 0.3 milliarcsecond in the visible. The 
addition of a Wollaston prism allows the study of polariza- 
tion across velocity dispersed systems, such as the circum- 
stellar disk of 7 Cas. The performance of the basic instru- 
ment is summarized by Mourard et al. (2011). 

Begun in 2008, the VEGA observations of 7 Cas summa- 
rized in this paper and as noted in Paper 2. As noted in the 
latter, half of them were contemporaneous with the XMM- 
Newton observations, and nost of these measurements were 
made over a variety of configurations. A limitation of our 
VEGA dataset is that the observations were clustered near 
the same minimum baseline (^30 m) in the South (S1S2) 
alignment. The W1W2 baseline was more than adequate 
to resolve the disk in the East- West direction on the sky. 
However, this mismatch in configurations limited our abil- 
ity to verify the presence of features in the South-North 
direction on the disk, even though departures from disk 
axisymmetry have been observed with earlier generations 
of equipment (e.g., Berio et al. 1999). Also, we were unable 
to get a complete set of phase closures at multiple times. 
This means that we could not determine time-dependent 
differences during the 2008-2010 interval. 

4.2. The Automated Photometric Telescope Observations 

Our optical continuum observations of 7 Cas were made in 
the standard Johnson V and B band niters on every avail- 
able photometric night (sometimes several times per night) 
using the 0.4 meter T3 Automated Photometric Telescope 
(APT) located at Fairborn Observatory in the Patagonia 
mountains in southern Arizona. Our monitoring of this star 
began in 1997 September and has continued every year 
since, with an observing season lasting from September into 
February and restarting again briefly in mid- June, weather 
conditions permitting. We lose about two months of the 
7 Cas observing season due to monsoon rains. 

Full details concerning our APT observational program 
on 7 Cas can be found in RSH and SHV. Briefly, the ob- 
servations are made by cycling the automated telescope 
between the program star and the nearby stars HR6210 
(F6V) and HD5395 (G8IIIb), which are used as compar- 
ison and check stars, respectively. The APT executes each 
observing sequence as follows: K,S,C,V,C,V,C,V,C,S,K, 



where K, C, and V, are the check, comparison, and pro- 
gram stars, respectively, and S is a sky observation. Each 
complete sequence is termed a "group observation" or sim- 
ply a group. All integration times were 10 s, and all mea- 
surements were made through a 3.8 mag neutral density 
filter to attenuate the counts and avoid large deadtime cor- 
rections. The program star is thus bracketed by the com- 
parison star on both sides three times during each group 
measurement at a cadence of 8 minutes. Group mean V — C 
and K — C differential magnitudes are computed and stan- 
dardized with nightly extinction and yearly mean transfor- 
mation coefficients determined from nightly observations of 
a network of standard stars. Typical rms errors for a sin- 
gle observation, as measured from pairs of constant stars, 
are ±0.0030 — 0.0035 mag. The rms scatter of the seasonal 
mean K — C differential magnitudes, inferred from differ- 
ences of the 40th and 60th percentile of fluxes during a sea- 
son (a reliable internal metric) is ~ ±0.003 mag., indicating 
good long-term stability. By comparison, the corresponding 
figure during the 2010 season was twice as large, ±0.007 
mag. This was also near the value of the formal rms. These 
larger values are largely due to the quasi-secular brighten- 
ing during this time. 

4.3. Optical wavelength spectra ofj Cas 

To support this program we made use of 306 echelle spec- 
tra of 7 Cas taken in a long-term monitoring program with 
the spectrograph fiber-fed from the Cassegrain focus of 
the 1.06-meter telescope at Ritter Observatory (adminis- 
tered by the University of Toledo) in 2009 December and 
the second half of 2010. The spectra were obtained using 
a new Spectral Instruments, Series #600 camera config- 
ured in an echelle format. This resulted in a resolution 
of 27,000. We used three echelle orders to extract the Ha 
spectrum and background, extending over the range 6285- 
6834 A. Exposures lasting a few minutes gave spectra with 
a signal to noise ratio of «100. 

Our dataset includes Ha spectra by E. Pollmann who 
used a cooled CCD camera with 9x9 /im pixels attached 
to the Cassegrain focus of his 35 cm Schmidt telescope in 
Leverkusen. A dispersion of 22.2 A mm made for a spec- 
tral resolution of 16,400 over the range 6520-6700 A. In 
Hamburg R. Biicke used a 20 cm Dobson Telescope fiber 
coupled to a spectrograph contain a cooled CCD camera 
with 18x200 (an pixels. This system provides a dispersion 
of 0.56 A pixel -1 and a spectral resolution of 3,400 extend- 
ing over a range 5900-6700 A. These instruments have been 
used in stable configurations for several years. 

4.4. XMM-Newton Observations 

We conducted four XMM-Newton observations under the 
auspices of Guest Investigator cycle A09. These were car- 
ried out on 2010 July 7 (ObsID 0651670201; referred to 
here as OBS1), July 24 (0651670301; OBS2), August 2 
(0651670401; OBS3), and August 24 (0651670501; OBS4), 
beginning on MJD dates 55384.30, 55401.15, 55410.77, and 
55428.09, respectively; XMM-Newton prefers the MJD time 
system. The observations OBS1 to OBS4 had total effec- 
tive exposure times of 17.5, 15.7, 17.5 and 22.4 ks for 
the pn detector, respectively, 0.2 ks longer times for the 
MOS1 and MOS2 cameras, and 0.4 ks for the RGS cameras. 
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High energy-sensitive data were taken through the thick 
filter of the European Photon Imaging Camera (EPIC) 
camera operating in small window mode. To avoid high 
noise contributions from soft energies, our timing analy- 
ses were restricted to EPIC data in the 0.8-10 keV energy 
band. Our spectral analysis made use of the complemen- 
tary information of the medium resolution EPIC spectra 
between 0.2 and 12 keV and high resolution RGS spec- 
tra from 0.3 and 2 keV. Data reduction, spectral analy- 
sis and part of the timing analysis were carried out using 
SAS vlO.0.0, XSPEC, and Xronos package, respectively. See 
http://heasarc.nasa.gov/xanadu/spec/ for a description of 
these packages. 

Altogether we used five XMM-Newton cameras in our 
observations. Of these the EPIC detectors (pn, MOS1, and 
MOS2) were especially important for the evaluation of the 
hot plasma component. The RGS1 and RGS2 cameras were 
more important for assessing the warm and cool plasma 
components defined below, and exclusively in determining 
spectral line attributes for lines above 10 A. The optical 
monitor was turned off during our observation. 

5. Results of optical investigations 

5.1. Fitting of CHARA visibilities of the 7 Cas system 

In this section we summarize the interferometric results ob- 
tained in our collaborative program. 

We were able to resolve the disk diameter in both ri- 
band and R-band using LITpro model programs (for de- 
tails see http://www.jmmc.fr/litpro_page.htm). We fit the 
MIRC data with two three-component models, that is with 
a uniform brightness star disk, a Be disk and a uniform 
"background" comprising 12-14 % of the total flux. Such 
a component is not the norm for Be-disk descriptions. 
Physically it could be due to the presence of previously 
ejected matter from the star that is now extended beyond 
the disk. In the first model the Be disk is assumed to have 
a Gaussian profile with distance from the center of the star. 
In the second model the disk model is an elliptical ring of 
equal surface brightness, with a brightness cutoff at the in- 
ner and outer edges. The motivation for the ring model was 
to determine whether there is a star-disk gap that might 
arise from solid body rotation of clouds within the radius. 
However, we were unable to get a satisfactory fit to the 
VEGA data using the ring model. Table 1 shows the re- 
sults. The angular values in milliarcsec here are converted 
to linear dimensions and represented in italics, assuming 
the distance and value of R* given in SJTJ 

Using these models with the MIRC data, we resolved 
a flattened circumstellar disk with elongation ratios of 
0.82/0.62 = 1.33 (Gaussian) and 0.85/0.61 = 1.39 (ring). 
The VEGA R-band data led to an elongation of 0.76/0.56 
= 1.36 The MIRC values give an inclination i«42° and 44°, 
respectively. The VEGA R-band results are essentially the 
same, iss42°±3°. Likewise, the position angle derived from 
the VEGA is 19°±5° and thus consistent with the MIRC 
value. For the MIRC solution the ring radius extends from 
1.4 R* to 2.0 R*, but the error bars are large. 

As Table 1 shows, we were unable to get a satisfactory 
solution for the VEGA R-band data with a ring model. 

We also used the VEGA system to measure fluxes in 
the Ha line both spatially and kinematically. In this anal- 
ysis we adopted the major axis and inclination angle from 



the R band analysis outlined above. Our best fitting ma- 
jor axis parameter in Ha flux is 4.4±0.4 mas. The position 
angle was measured as 20±15°. Assuming the previous pa- 
rameters, we find a Be disk radius of 10±1R*. We point 
out that at 10 R* the VEGA disk radius is nearly half the 
Be star's Roche lobe radius (of some 21.4 R*; Gies et al. 
2007). We find from Okazaki & Negueruela (2001) that the 
most likely 3:1 resonance truncation radius of the disk due 
to the companion in a circular orbit is about 0.4 of the 
orbital separation, or « § of the Roche lobe radius. Our 
VEGA radius does not yet reach this extent, but one may 
hope more sensitive future observations, particularly in the 
far-IR, that could test for the presence of a truncation edge. 

Our kinematic observations for Ha were designed to 
determine the exponent value (3 in the velocity distance 
relation in the disk, = (r/R*)' 3 . We obtained a value 
of /3 = 0.5±0.1, or exactly the Keplerian rotation result. 
Solid body rotation, for which j3 = -1, is therefore (un- 
surprisingly) ruled out for the whole disk. The confluence 
of Keplerian rotation with corotating clouds over the star's 
surface suggests a complicated geometry for the CS matter. 

5.2. Optical photometry and spectroscopy 

5.2.1. Combining the APT light curve and Ha measures 

FigurefTJi plots the seasonal mean Johnson V magnitudes 
and mean (B — V) color indices from the T3 APT for the 
13 observing seasons from 1998-20100 

Notice a series of long-term erratic undulations over 
timescales of 2-3 years. This variability pattern probably 
results from structural changes in the disk, but the details of 
this are unclear. We can shed some light on this by coplot- 
ting the Ha measures in Figured)}. Although the light curve 
shows highly statistical fluctuations from year to year, over 
the entire timespan the slopes of both Ha equivalent width 
and F-band light curves are positive - that is, the activ- 
ity as defined by Ha emission, magnitude brightening, and 
B — V color have increased more or less together. Thus 
the Ha emission from the CS disk correlates with optical 
continuum brightening. 

Inspection of the two curves in Fig. [1] shows that second 
order deviations in their behaviors exist. These include a 
delay in the decrease in the Ha emission index over the 
Reduced Heliocentric Julian Date (RHJD) range 51040- 
51080, absences of rises in the index at RHJD 52000-52400 
and at 54000-54200, and absence of a rapid upturn at the 
end of the data sequence. Such deviations can arise because 
of differences in the optical thickness in the V^-band con- 
tinuum and the Ha line core. The optical depth in the Ha 
line is many times more optically thick than the y-band 
continuum (see e.g., Carciofi 2011) and therefore the point 
at which Tdisk = 1 lies further from the star. Consequently, 
the Ha emphasizes disk changes at greater radius than the 
V continuum does. 



3 We define the "1998" season as the interval from 1998 June- 
1999 February. Note that we omit the 1997 observing season 
because a neutral-density filter change in 1998 resulted in new 
magnitude zeropoints from 1998 onwards. Analysis of the vari- 
ability in 1997 was published by SHV. Errors in the seasonal 
means are very small and are due mainly to incomplete averag- 
ing of the long cycles and/or slopes in the magnitude. 
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Table 1. Relevant star and Be disk fitting parameters from LBOI 



MIRC DATA: IR H-band 



PARAMETER 


GAUSSIAN DISK MODEL 


RING DISK MODEL 


Star Diameter (mas) 


0.44±0.03 


0.44±0.03 


Inclination angle (i) 


42±1° 


44 ±1° 


(thin disk or ring) 






Minor disk axis (mas) 


0.62±0.08 




(Linear, in 2Rq ) 


1.4 




Major disk axis (mas) 


0.82±0.08 




(Linear, in 2Rq ) 


1.86 




Minor ring inner diam. (mas) 


n/a 


0.60±0.33 


(major axis) 






Outer ring diameter (mas) 


n/a 


0.85±0.45 


along the major axis 




Flattening ratio 


1.33±0.08 


1.39±0.08 


Position Angle 


12±9 


12±9 


(E from N in degrees) 






Star flux (%) 


35±1 


45±6 


Disk/ring flux (%) 


53±2 


41±6 


Background flux (%) 


12±1 


14±1 


x 2 


3.96 


3.95 


VEGA: near-Ha "R"-band 


PARAMETER 


GAUSSIAN DISK MODEL 


RING DISK MODEL 


Star diameter (mas) 


0.44±0.03 


(no fit) 


Inclination angle (i) 


42±3° 


n/a 


(thin disk) 




Minor disk axis (mas) 


0.56±0.05 


n/a 


Major disk axis (mas) 


0.76±0.05 


n/a 


Major axis (in 2Rq) 


1.72 


Flattening ratio 


1.36±0.08 


n/a 


Position Angle 


19±5 


n/a 


(E from N in degrees) 




Star flux (%) 


55±2 


n/a 


Disk flux (&) 


45±2 


n/a 


x 2 


1.95 


large 



5.2.2. A revised radial velocity solution for the binary orbit 

We measured the radial velocities for the H a spectra from 
Ritter Observatory by computing bisectors of the emission 
line components. For this purpose we used the method of 
Shatter, Szkody, & Thorstenson (1986), which has special 
application to measurements of shifts of asymmetric line 
profiles. The bisector position is determined by creating 
a function composed of a negative- valued Gaussian (offset 
by AV = —300 km s _1 ) and a positive-valued Gaussian 
(offset by AV = +300 km s^ 1 ); both Gaussians have a 
width FHWM = 85 km s _1 . This function is then cross- 
correlated with the observed emission lines (whose base- 
line is reset to zero), and the zero-crossing of the result- 
ing cross-correlation function yields the bisector position. 
We measured 198 velocities from Ritter Observatory that 
span from 2000 January 28 to 2010 November 11. Also, 
one of us (RB) obtained 103 new measurements over the 
interval 2006 September 20 through 2011 April 22 from 
the Hamburg Dobson Telescope system. These were mea- 
sured by fitting a Gaussian function to the profile. Finally, 
we combined the Ritter and RB velocities. The weights of 
the Ritter data were determined from the mean continuum 
signal-to-noise ratio in the spectra. Because the RB spectra 



Table 2. Ritter & Dobson/RB radial velocities (stub) 



RHJD 


ARV 


Error 


O-C 


Phase 


Source 


51571.510 


1.63 


0.92 


-2.16 


0.00 


Ritter 


51586.504 


3.04 


1.32 


-0.30 


0.08 


Ritter 


51608.527 


0.38 


1.31 


-1.09 


0.19 


Ritter 


51609.516 


-0.22 


1.69 


-1.59 


0.19 


Ritter 


51777.805 


3.80 


5.81 


0.04 


0.02 


Ritter 


51794.820 


2.87 


1.64 


-0.17 


0.10 


Ritter 


51800.800 


0.31 


1.19 


-2.26 


0.13 


Ritter 



were measured from lower resolution spectra, we assigned 
them a higher velocity error of ±3.16 kms -1 . The RB ori- 
gin can be found from this assignment in Table 2. To bring 
the Ritter and DB velocities to a common mean of zero, we 
added 5.73 kms -1 and 9.95 kms -1 to the Ritter and RB 
velocities, respectively. The results are shown in Table 2 as 
a stub (dates in Reduced Heliocentric Julian Dates). The 
full data list of 301 radial velocities is given in the elec- 
tronic version of this paper. Column 3 is the measurement 
error in kms -1 ; column 4 is the "O - C" (Observed minus 
Computed) differences; column 5 the orbital binary phase, 
and column 6 is the RV source (Ritter or RB). 
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Table 3. Orbital Solution for 7 Cas. 



Parameter 


Units 


1 


2 (adopted) 


K 


(km s" 1 ) 


3.80 ± 0.08 


3.79 ± 0.08 


e 


0.03 ± 0.02 




Ul 


(degrees) 


198 ± 29 




T 


(RHJD) 


55549.48 ± 16.32 




JilVmax 


(RHJD) 


55436.47 ± 0.64 


55438.1045 ± 1.35 


P 


(days) 


203.590 ±0.21 


203.555 ±0.21 


f(M) 


(M e ) 


(0.116 ±0.008) x 10~ 2 


(0.115 ±0.007) x 10~ 2 


ai sin i 


(Re) 


15.28 ±0.33 


15.23 ±0.33 


rms 


(km s" 1 ) 


1.03 


1.03 




1998 2000 2002 2004 2006 2008 2010 2012 
Calendar Date 

Fig. 1. (a): APT V (my) magnitudes and (B — V) color of 
7 Cas over the 1998-2010 seasons of coverage; annual aver- 
ages are shown. Error bars on these magnitudes (marked as 
* and X, respectively) are much smaller than the symbols, 
as each point represents a few hundred individual observa- 
tions. Horizontal bars denote the coverage in time for each 
season's average. Small points indicate each my value mea- 
sured in the 2010-2011 season. Annotations at the bottom 
give the column densities n^ in units of 10 22 cm~ 2 respon- 
sible for the attenuations of soft X-rays in high dispersion 
X-ray observations of 2001, 2004, and 2010. (b): The Ha 
emission equivalent widths (EW) in Angstroms obtained 
by Pollmann and Biicke (star and square symbols, respec- 
tively over the same time as panel (a)). A displacement of 
-1.0 A was required to bring the the mean Biicke EW to 
the Pollmann mean. The vertical dotted lines correspond 
to times of the first and last XMM-Newton observations 
in 2010. Note the rapid Ha emission increase from 2002 
on, and also a similar spike in 2007 as observed in 2010. 
Both events coincide with my brightenings when panel a is 
plotted in high time resolution. 



In order to construct a new orbital solution, we first ex- 
amined these new data and found evidence for long-term 
trends introduced by including pre-2000 radial velocities. 
We did not include these early data because the methods 
for pre-whitening could introduce new artifacts in the anal- 
ysis of the orbital elements. We made an error-weighted fit 
to the velocities to solve for orbital elements using the pro- 
gram of Morbey & Brosterhus (1974). We solved for the 
orbit by allowing the eccentricity e as well as other orbital 
parameters to float freely. This solution is represented as 
Orbit 1 in Table 3. Discovering that this eccentricity is at 



o 
o 

CD 

> 



"O 

CO 




0.5 

Orbital Phase 

Fig. 2. The adopted radial velocity solution from our 
Orbit 2 model (see Table 3). Solid dots and open trian- 
gles are measures obtained from Ritter Observatory and 
RB's Hamburg system, respectively. The weights to the si- 
nusoidal solution are lower for the open symbols. Table 2 
gives errors for all observations. 



most no more than marginally significant (e.g., from testing 
the probability index p\ suggested by Lucy 2005), we ulti- 
mately adopted the circular solution with the same period. 
This is given Orbit 2 in our table. We depict the folded data 
to the Orbit 2 solution as a 203.555 day sinusoid in Figured! 
Phase zero for Orbit 2 is referred to the maximum velocity. 
For Orbit 1 we adopt the time of periastron passage. 

Our derived orbital parameters for 7 Cas differ slightly 
from the discovery analyses of Harmanec et al. (2001) and 
Miroshnichenko et al. (2002). The first study found e = 
0.26 and P or f, = 203.59 days, and the second e = 0.0 
and P rb — 205.50 days. Essentially, our new period solu- 
tion agrees with Harmanec et al. but our eccentricity value 
agrees with Miroshnichenko et al. Shortly before submission 
of our results for publication, we received a preprint from 
J. Nemravova and collaborators describing a radial velocity 
study covering nearly 17 years of the star's velocity history 
(Nemravova et al. 2012). Their results and ephemeris are 
substantially the same as our own. Their estimate of the 
period, 203.52±0.08 days, is well within the formal errors 
of our two analyses. The same comment applies to the er- 
rors in To and phases. Therefore, we believe a synthesis of 
our works leads to a definitive orbital parameters. 

From the parameters in Table 3, if we assume a mass of 
15M0 for the primary and an inclination of 45° and assume 
that the binary orbit and Be disks are coplanar, we may 
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adopt the mass function of 1.15x 10 _3 Mq to find a most 
probable secondary mass of 0.8 Mq. Assuming reasonable 
errors for the primary mass (±1 Mq), an inclination of 45° 
(±6°), and Ki (±0.1 kms -1 ), we estimate a generous error 
of ±0.4 Mq for the secondary. This estimate supports the 
possibility that the secondary could be a late type dwarf, 
white dwarf, or a sdO star, but probably not a neutron 
star. This justifies the omission of this star from adding an 
additional optical stellar component in our LBOI models. 

If we adopt the Orbit 2 determination in Table 3, our 
orbital phases corresponding to the XMM-Newton observa- 
tions discussed below are (f> or b — 0.74, 0.82, 0.87, and 0.96 
(where again 4> or i, — 0.0 is defined relative to the epoch of 
maximum RV). The Nemravova et al. RV orbital ephemeris 
gives phases for these observations of 0.73, 0.81, 0.86, and 
0.95, i.e., within 0.01 of our values. As for our solution, 
Nemravova preferred a circular solution to even small val- 
ues of the eccentricity. In considering the relative weighting 
of these two RV studies, we can take the errors at face value 
and double weight the Nemravova et al. solution relative to 
our own. This gives a period of 203.53±0.10 days. 

To summarize, the first of the 2010 XMM-Newton ob- 
servations was conducted just as the Be star was in infe- 
rior conjunction with respect to its secondary, similar to 
the phase of the 2001 Chandra observation. The last of the 
four observations coincided with the Be star essentially at 
quadrature phase. With this emphemeris one finds that the 
phase of the 2004 XMM-Newton observation was 0.27. 



6. Analysis of XMM-Newton Observations 

Herein we discuss four high-quality XMM-Newton expo- 
sures of 7Cas in 2010, dubbed OBS1-4 above. The star 
was previously observed in high resolution on two occa- 
sions: in August, 2001, by the High Transmission Energy 
Grating of Chandra (HEG/MEG) and in February, 2004, 
by XMM (RGS) (Smith et al. 2004, Lopes de Oliveira et al. 
2010). We can combine the photometric and spectroscopic 
results to get a sense of epochal differences of the conditions 
responsible for X-ray emission. 



6.1. XMM-Newton light curve variations 

We began the analysis by extracting light curves binned 
to 20 seconds for hard- (4.5-10 keV), soft- (0.3-1 keV), and 
the total flux over 0.3-10 keV for each of the four obser- 
vations. The median fluxes of three of the light curves ex- 
tracted from the pn detector (0.3-10 keV) are about the 
same, about 25.0 counts s _1 . The mean flux of the third 
observation at MJD 55410 is 47% higher than the mean of 
the other three, among which the difference is only ±7%. 
Given the assumed distance of 168 pc, this means 
that the L^'s for the observations lie in the range 
8— 13xl0 32 ergs -1 . This range is somewhat higher than 
the range 3-9 xlO 32 ergs -1 from RSH using the RXTE 
several years earlier. 

The light curves and hard/soft color ratios for OBS2 
and OBS4 are shown as examples of all four observations 
in Figure|31 Note that the light curve exhibits the usual 
pattern of almost uninterrupted flaring. 

As with previous light curve studies, we find also that 
soft- and hard-energy light curves constructed track one 
another very well, giving only occasional significant depar- 



4.5-10 keV cur\ 



XMM OBS 2 



^4.5-10/0.3-1 ratio 



—4.5-10 keV curve 
25 '-. I 



g ^4.5-10/0.3-1 ratio Ratio 



Time after start (ksec) 

Fig. 3. Two stretches of XMM-Newton light curves and 
color ratios for OBS2 and OBS4. The light curves are 
binned to 20 seconds and smoothed over two points. The 
rms errors are indicated taken from point to point fluctua- 
tions, with the same smoothing. The color ratios have been 
scaled by 5 x for visual clarity. Occasional flare aggregates 
are obvious in the light curves. The lower plots in the two 
panels disclose the color variations only over timescales of 
J>0.6 kiloseconds. 



tures from a mean color ratio. There is in general a remark- 
able absence of color hardening during strong flares or brief 
near cessations in emission. In fact, most of the largest color 
fluctuations occur when no strong flux excursions (flares or 
flux lulls) occur. Yet we do find occasional departures from 
a uniform color. For example, we found a few sinuous vari- 
ations in color curves in each of the four observations. An 
example is depicted for stretches for OBS2 and OBS4 - see 
Figure[31 In OBS4 sudden hardenings occur over 0.5-1 ksec 
and are followed by a slow return to the mean color over a 
few ksec. Sometimes hard flux changes are followed by soft 
fluxes, but not quite to the same degree (e.g., the "decays" 
in Fig.[3jD). At other times they occur when the hard com- 
ponent is nearly constant or undergoes apparently stochas- 
tic fluctuations (Fig. (3^.). In almost all these cases the soft 
component changes more slowly, on timescales of 100-1000 
seconds, in the fashion shown in these plots. 

The power spectrum, summed over our four observa- 
tions, is exhibited in Figure[4] It shows the beginning of a 
logarithmic slope just exceeding -1 for > 30s flare-like ex- 
cursions. At about 0.005 Hz it breaks to a shallower slope. 
This is consistent with the break at 1-3 xl0~ 3 Hz reported 
by Robinson & Smith (2000). We characterize the occa- 
sional variations on timescales longer than flares as "un- 
dulations." SRH and RS have shown that these features 
are most likely due to the passage of translucent clouds 
forced into corotation over the Be star. They suggested that 
these are responsible for the continued rise in power above 
the break in the slope of the power spectrum. Looking to 
shorter timescales, we performed an autocorrelation anal- 
ysis on the integrated light curves as well as light curves 
separated into binned soft and hard fluxes. However, we 
found that none of them exhibited quasi-periods. Also, we 
found that to a precision of 1-2 sec we could find no evi- 
dence of time delays between the hard and soft flux time 
series. 
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Fig. 4. Sum of power spectra of our four /XMM/EPIC 
pn observations over 0.8-10 keV of 7Cas. The binning in 
time over 10 seconds prevents the spectrum from being fol- 
lowed out to its white noise component. Regression lines 
are drawn with slopes noted through two frequency ranges. 
A break in the slopes is noted by a square symbol at a 
frequency of 0.005 Hz. 



6.2. Spectral analysis: primary parameters 

Description of XMM-Newton spectra 

The 2010 XMM-Newton spectra appear similar to those 
obtained by Chandra in 2001 (Smith et al. 2004) and 
the XMM-Newton in 2004 (Lopes de Oliveira et al. 2010, 
"L10"). The spectrum is "hard," i.e., most continuum flux 
is emitted at energies of a few Angstroms or lower. The long 
wavelength continuum is attenuated much like the Chandra 
spectrum and in contrast to the unattenuated 2004 XMM- 
Newton spectrum. The visible emission lines in the long 
wavelengths of the RGS spectral detectors are the Lya- 
like lines of even-nucleons MgXII, NeX, and OVIII, and 
the i and / components of the so-called fir complex of the 
ground- level helium-like ions Ne IX and O VII. The Lya-likc 
emission lines of the Fe K shell ions Fe XXVI and Fe XXV 
are likewise visible in the short wavelength spectra of the 
EPIC/ MOS and pn detectors. An additional emission fea- 
ture at wl.94 A arises from atoms of less excited Fe ions, 
the so-called Fe K fluorescence feature. Similarly, in three 
of the four RGS camera spectra we found faint correspond- 
ing features due to SiK ion fluorescence, as previously re- 
ported by Smith et al. (2004, "S04"). In addition, the RGS 
spectra exhibit a number of emission lines of Fe XVII and 
Fe XVIII (Fe L shell) ions. The quantitative analyses of 
our line strengths and broadenings were performed inde- 
pendently for spectra recorded by the pn and MOS detec- 
tors at short wavelengths and for the RGS\ and RGS2 
detectors at long wavelengths because the instrumental re- 
sponses and calibrations for these cameras are different. 
Following XMM-Newton mission recommendations we have 
performed quantitative analyses with XSPEC using all in- 
dividual EPIC and RGS spectra simultaneously. 

Using XSPEC to determine multi-temperature models 

We used the same suite of ftools and XANADU '/XSPEC 
tools to fit our XMM-Newton spectra that L10 used - 
mainly mekal, vmeka, and cevmkl. We also used the XSPEC 
model bvapec to determine measures of line widths and 



certain elemental abundances. All of these models assume 
one or more optically thin and thermal plasma components 
that are also in coronal ionization equilibrium. The mekal 
model families are then parameterized by a series of mono- 
thermal plasma components and their normalizations (pro- 
portional to Emission Measures, "EMs"). The normaliza- 
tions in the various cameras were independent of one an- 
other and all were kept free during the fit. This procedure 
was adopted to assure that each individual normalization 
takes into account the difference in calibration of the cam- 
eras. Ultimately we adopted the mean value of the normal- 
ization of each camera and for each plasma component to 
compute an emission measure from its recorded spectrum. 
However, in doing so we neglected the weak and uncer- 
tain high energy response of the RGS cameras, which can 
otherwise bias the result. A Gaussian emission feature was 
added to the model to parameterize the visible Fe fluores- 
cence line at 1.9 A. We also included necessary multiplica- 
tive photoelectric absorption terms, parameterized by the 
term "nj^ ," which correspond to the equivalent hydrogen 
columns. These parameters are combined in a simultaneous 
solution to fit the continuum and emission lines. In general, 
a number of iterations are required for convergence. We ex- 
perimented with different starting temperatures to avoid 
converging to false solutions. 

As in L10 we tried to fit the 7Cas XMM-Newton 
spectrum with a continuous Differential Emission Measure 
(DEM) defined as a cooling flow model (e.g., cevmklf^. 
While this is sometimes considered the preferred X-ray 
emission description for accreting white dwarfs, it does not 
work for our spectra. For example, for adiabatic flows and 
Tmax allowed to float to 28 keV, the continuum fit is poor, 
giving a reduced x 2 — 1-69. For a model where the a index 
is permitted to float to its optimum value of 1.31, the fit is 
still poor, x 2 = 1-52. 

Next we ran multiple-component optically thin thermal 
models with mekal for the five XMM-Newton spectroscopic 
cameras. The final results are shown in Table 4. These mod- 
cling results were similar to the results of S04 and L10 in 
that they led us to an optically thin model with thermal 
plasma components and two absorption columns. The so- 
lution for the cool plasma component denoted by kT2 is an 
increase of 0.8 keV from the L10 model - a statistically sig- 
nificant result due to an actual spectral change according to 
our return to the L10 spectra with our mekal models. The 
addition of successive temperature components reduced the 
mean x 2 from 1.70 to 1.48, to 1.28, and to (four compo- 
nents) 1.20. With each added component the improvement 
in fit over certain wavelength regions was marked. In par- 
ticular, mekal models with fewer than four components also 
caused perceptible spurious undulating departures from the 
observed continua that were minimized with the inclusion 
of the fourth component. The addition of further compo- 
nents reduced the mean x 2 by negligible amounts and were 
not significant according to statistical / tests performed in 
XSPEC. However, like L10 we were not able to fit both the 
pn/MOS and RGS spectra well unless we added a second 
column to attenuate the soft energies. We may parameter- 
ize our best mekal model with the equation: 



4 The differential emission measure for cooling flows is given 
by dEM = [(T/T max ) a ](dT/T max ). For adiabatic flows a = 1. 
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Fig. 5. Color-coded "unfolded" model spectra in the range 10-30 A computed for RGS1 and RGS2 for the fitted model 
of OBS3. The green, blue, purple, and red denote the normalized emission spectra of plasma components fcTi, fcT2, fcT3, 
and fcT4, each with the single values given in Table 4. The black line is the observed spectrum. The model from the 
RGS2 detector was used for all wavelengths except for the 20-24 A gap, and for this gap the RGS1 model was used. This 
accounts for small discontinuities in the emissions at these wavelengths for the fcT2 and fcT3 plasma components. The 
/CT3 and fcT4 emission are almost as well described by a broad DEM centered at a high energy. 

Table 4. Best parameters from X-ray Model of the form: (Ti +T2+T3 + T4) * riH a + (T4 + gauss) * n.H b - 



Parameter 

n Ha (10^ cm-' 2 ) 

fcTi (keV) 

EMi (10 55 cm" 3 ) 

kT 2 (keV) 

EM 2 (10 55 cm" 3 ) 

fcT 3 (keV) 

EM 3 (10 55 cm" 3 ) 

fcT 4 (keV) 

EM 4( , (10 55 cm" 

EM 4 , . (10 55 cm" 
b 

f T (erg cm" s" ) 
n Hb (10 22 cm" 2 ) 
redshift 

Xg/d.o.f. 



OBS1 



OBS2 



OBS3 



OBS4 



n 100+u.uul 

u ' iJO -0.001 


u - lzy: -o.ooi 


090 +u - uui 


198+U.UUl 
U - i28 -0.001 


n iir+0.002 

U.llO_ 002 


n -,-, n +0.002 
U.11U_ 002 




111+0 002 
u - 111 -0.002 


0.14 


0.10 


0.21 


0.12 


0.60_°;°_ 


fi4+ 02 


0.63ir o i 


o-62±S:Si 


0.33 


0.27 


0.37 


0.32 


O /?q + 0.22 

°- UJ -0.19 


4 ™ +0 ' 22 


4 20+ ' 14 


3.43l^ 3 


0.18 


0.29 


0.45 


0.22 


14 66 +0 - 31 


1 4 9Q+ - 25 
14._o_ 2 2 


1 K 71 +0-22 
10. 1 l_o.23 


1Q ci +0.18 
ro.or_0.i6 



3.0 
0.7 

2.63xl0" 10 

o.oois 

1.15/3189 



2.6 
1.0 



2.49x10 
57.0±42 
0.0005 
1.18/2877 
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3.5 
0.8 

3.89xl0" 10 

59.6^26 
0.0007 
1.18/3663 



2.6 
0.7 



1.79x10 
73.7± 3 ; 2 
0.0003 
1.20/3441 
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The following elemental abundances were fixed to nonsolar values: ZFeKa = 0.18xZ_F e ,Q, Zn = 2.33xZat,q, ZNe = 1.8xZjv e ,0. 
Nonzero redshift values are not considered statistically significant because they are of the order of known seasonally-induced 
shifts in the instrument (e.g., Pollack 2010). Entries marked "fT" are unabsorbed fluxes in the energy range 0.2-12 keV. 



phabs a * (mekali + raekali + mekal^ + mekal4 : ) + 

+ phabsb * (mekal^ + gaussian) , (1) 

where the mekal subscripts 1-4 parameterize the tempera- 
tures fcTi-fcT4 of the four thermal plasma components of 
a 4-T model. The subscripts in phabs express the two in- 
dependent column densities, rin a and nu b . As we will see, 
these values can be very different. We stipulate that it is 
possible in principle that the fcT2- 4 components could be 
represented in the nn b term as well, but their fluxes are so 
much weaker that we cannot solve for them as additional 
independent variables. In Figure[S]we exhibit the unfolded 
model solution, including the normalizations, as an exam- 
ple for the range of 10-30 A in OBS3. This gives a good 
feel for the relative contributions among the various plasma 
components for both the line and continuum spectra. 

Figure[5] shows that even in the soft X-ray region the hot 
kT4 component dominates the flux. In t )6.1l we noted that 
the hard/soft ratios can change over timescales of 100-1000 



seconds. Despite the dominance of the hot emission com- 
ponent in the soft regime, the cool kTi plasma can some- 
times make its presence felt independently of the hot one. 
The most likely means for producing soft X-ray variability 
is for the line of sight column of matter to change to the 
hot emission sites behind it. If one assumes that the hot 
component fluxes decreased by free cooling, i.e., propor- 
tional to the reciprocal of the density, this timescale gives 
an electron density of ~T0 n cm" 3 . 

Whereas components kT2 and fcTi are critical to the 
fits of lines of the Fe L-shell, N, and C ions, the /CT3 com- 
ponent is the least stable in our solutions. One reason for 
including fcT3 is that it easily produces the FeXXV Lya 
line and contributes to part of the continuum. However be- 
cause the continuum flux spectrum of our fcT3 runs paral- 
lel to the fcT4 in a logarithmic plot (Fig. [5]), it is possible 
that this component is an artifact of what is actually a 
continuous Differential Emission measure (DEM) parame- 
terized by AT4 but extending down to a few keV. To test 
this possibility we ran an XSPEC add-on package written 
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by K. Arnaud called vgadem, which is a DEM response 
shaped like a Gaussian, with a peak emissivity at 12.25 keV 
and a sigma of 5.47 keV. The values of the /cTi and fcT2 
were essentially identical with those determined with the 
mekal model. The resulting model replaces the monoener- 
getic emissivity terms, fcT3 + kT^, in the parenthesis term 
of equation (1) and kT^ in the second term with one having 
a broad DEM. For our test case (OBS1) this model gave an 
acceptable \ 2 °f 1-19 f° r the global fit. However, we found 
this solution to be unstable in convergence, and it slightly 
underpredicts the strength of the FeXXVI line. In all, it 
does not perform as well as the mekal 4-T model over all 
wavelengths. We had a similar result running vgadem with 
the 2004 spectrum and abandoned this model. From these 
trials we could not corroborate our hypothesis that the kT^ 
plasma is an unnecessary product of our reference model. 

We pursued this investigation by comparing our mekal 
4-T and 3-T models and excising the wavelengths around 
the FeXXV line. We found that the inclusion of the 
apart from these wavelengths causes the 4-T model to fit 
the continuum slightly better (x 2 from 1.19 to 1.15). Thus, 
the kT^ in mekal seems needed for an ideal fit to the con- 
tinuum too. However while introduction of the kT^ im- 
proves the fit for FeXXV in our reference model, the fit 
for FeXXVI is actually made slightly worse than for the 
reference model. To conclude, although we cannot provide 
an objective reason for preferring the 3-T vgadem DEM 
model over the reference model, the mekal model with fcT3 
* 4 keV actually appears not to be ideal either. We suspect 
a flatter DEM with cutoffs at low and high energies could 
improve these fitting deficiencies. In fact, there are phys- 
ical reasons for believing that the high energy source for 
the X-ray emissions of this star is not single- valued in en- 
ergy even as there are reasons for suspecting that the "A:T3 
component" is in fact not independent of the kT^ one. 

We give the parameters found in the mekal solutions 
represented by equation (1) in Table 4. The solution also 
stipulates the anomalous abundances found. Some 80% of 
the flux is emitted by a hot component we designate kT^. 
This percentage increases to about 87% for the case that 
/CT3 represents an extension of &T4. 

The &T4 component dominates emission at all wave- 
lengths, so much so that it serves to weaken color ratio 
changes represented in Figure[3] when the kT^ component 
alone does not vary. This emission also dilutes the equiva- 
lent widths (EWs) of some lines produced in cooler plasmas. 
Figure[5] shows the unfolded model RGS spectra for OBS4 
from the 4-temperature model described in Table 4. We re- 
mind the reader that the &T4 flux is composed of a flaring 
and a dominant basal component. 

6.3. Analysis of secondary parameters 

Column densities, abundances, and line broadening: 
We list additional details emerging from our analysis: 

— the nn b values are the largest documented to date and 
are a few times larger than the 1 x 10 23 cm -2 found by 
S04. The values change by «2x among the four spectra. 

— the ratios of volumes associated with the "lightly ab- 
sorbed" (riff a ) and "heavily absorbed" (n,H b ) columns 
(given by the respective EMs in Table 4 and in tables in 
S04 and L10) remain the same during 2001, 2004, and 



2010. From Table 4 these ratios are 25±J 5 %, 23±3%, 
and 27±5%, respectively. 

— the bvapec task was used to determine abundances of 
elements from all the strong lines. Like several previous 
studies, the Fe abundance derived from the L shell ions 
is close to the solar value. The abundance determined 
from the two Fe K shell ion lines is 0.18x the solar 
value. This is statistically significantly higher than the 
even lower value of 0.12 found by L10. We used bvapec 
also to determine Ne and N abundances of 1.80 and 2.33 
times solar, respectively, from their Ly a line strengths. 
This is slightly more moderate than the 2.63±0.78x 
and 3.96±0.28x solar values that L10 found. We found 
no statistically significant changes in these abundances 
among the four 2010 observations. 

— the XSPEC model bvapec was used on the strongest line 
profiles in the RGS spectra to determine Gaussian line 
broadening. Our solutions were dominated by the strong 
O VIII line. We found FWHM values in the range 540- 
950 kms -1 . These values are marginally larger than the 
moderate values of 400 km s _1 and 478 km s _1 found 
by S04 and L10, respectively. 

— FeK and (barely) SiK fluorescence features are present. 
The EW of the Fe feature is -35mA to -50 mA. 
Although these features were present in earlier high res- 
olution spectra, their strengths in the 2010 spectra are 
the highest reported so far. 

Spectra derived from low and high flux time intervals: 

In order to understand how key fitting parameters 
change when 7 Cas's emission changes between low to high 
states (usually due to flaring), we split up the good time 
intervals according to low and high fluxes (relative to the 
means) and reconstructed spectra into low- and high- flux 
pairs. We discovered that the low and high flux times were 
strongly grouped together, creating split segments of times. 
The duration of these various segments was from 3-30 min- 
utes. This timescale is longer than durations of flares. 

The flux differences between these pairs were statisti- 
cally significant, but only for the kT^ plasma component. 
Except for OBS1, the high-flux spectrum was several per- 
cent hotter than the low-flux one. The largest difference 
in the solution for the /CT4 component occurred for OBS3. 
In this case the kTi values for the high and low flux val- 
ues were 16.08 keV and 14.59 keV, which, given the errors 
of ±0.59 keV, is statistically significant at more than two 
sigma. (We repeated this experiment for the 2004 XMM- 
Newton observation by splitting it into two equal time seg- 
ments and examining the low and high flux states in each 
half. We found no evidence for temperature differences for 
this observation.) Differences in the column densities be- 
tween the pairs were not significant. 

Even so, we found several statistically significant differ- 
ences among the high-/low-flux spectra among several emis- 
sion lines, most especially those arising from the Fe L shell 
ions, an emission feature at 16 A (the wavelength of O VIII 
Ly j3) and a neighboring NeX line, and also a blend of the 
NVII Lya and a nearby NVI line at 24.8 A. The general 
pattern is that these lines are stronger in the low-flux spec- 
tra. To shed light on these changes, we constructed unfolded 
spectra for each of the 8 high/low-flux spectra (two for each 
one of the four observations), according to their normal- 
izations. Examples of this unfolded low-flux spectrum for 
OBS2 and OBS3 are shown in Figure[5] The prominence 
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IS 15.5 16 16.5 17 

Wavelength (A) 

Fig. 6. Unfolded spectra for the spectral region shown in 
Fig-El f° r OBS3 (same color coding as in Fig. [5]). Most of 
the lines arise from the L-shell ions of Fe XVII & XVIII. 
Note that the strengths of most lines depends on the nor- 
malization of the kT2 emission component (blue) relative 
to the dominant kT4 component (red). 

of these lines was actually due to a diminished dominance 
of the fcT4 continuum, but not the lines, for the low spec- 
trum fluxes. In other words, the spectra of all the cooler 
components then "peek through" the dominant continuum 
emitted by kT^. The result is a highly dynamically chang- 
ing spectrum. The details of the low/high flux spectrum 
differences can be summarized in the following: 

— differences in strengths of most lines but OVIII, e.g., 
Fe L shell ions of FeXVII and Fe XVIII, NeX, and 
FeXVIII-blended Ly/3/Lya, especially in OBS2 and 
OBS3. The change is the smallest for OVIII Lya be- 
cause the kT2, fcT3, and kT^ plasma components con- 
tribute comparable amounts of EW. Therefore, changes 
among these components' normalizations, for example 
in the spectra extracted from low and high fluxes, do 
not affect their equivalent widths significantly. 

— analysis of unfolded spectra constructed from the high 
and low fluxes shows why these rapid variations in line 
strengths occur: the kT^ emission component in the low- 
flux spectrum decreases and allows emission from pri- 
marily the kT2 component to be more visible. An ex- 
ample is displayed in Figure[7p and modeled in Fig. [6] 
The change in riH b for low/high (OBS3) is of the same 
order as the deviations around the observational aver- 
ages. 

Other lines and density diagnostics: 

Similar to S04 and L10, where the fir (forbid- 
den/intercombinational/radiative) triplet complex is ob- 
served, in this case for NVI and OVII, we find that the 
ratio G = (i + f)/r ~ 1, and this means that the dominant 
process for formation of these triplet lines is consistent with 
the dominance of collisions in determining the relative level 
populations. Likewise, we find the ratio f/(i + r) ~ 0. This 
might be explained by the metastability of the upper level 
(ls2s 3 Si) of the forbidden transition. However, it is at least 
as likely that this level is radiatively depopulated by intense 
UV flux from a nearby hot star, in which case the utility 
of the forbidden component as a density diagnostic is vi- 
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Fig. 7. XMM-Newton/RGS spectra obtained from high 
(dashed line) and low (solid) flux good time intervals 
in the 15-17 A (FigJT]) region in which several density- 
sensitive FeXVII lines are present for OBS2 and OBS3. 
The generally greater strengths of these lines and the 
16.8A/17.0A ratio are discussed in the text. These spectra 
are smoothed over two 2 mA bins. 

tiated. Thus, as with the former studies, we cannot draw 
conclusions about density from these features. 

Other density-sensitive diagnostics are available among 
several Fe XVII lines at wavelengths accessible to both the 
Chandra/HTEG and it XMM-Newton/RGS wavelength re- 
gions. These have been studied empirically and theoreti- 
cally in detail by Mauche and collaborators (Mauche et al. 
2001, 2004) in the well observed Chandra/HTEG spectra of 
the mCV (polar) star AE Aqr. The lines in this source ap- 
pear to be emitted in plasma with thermal conditions sim- 
ilar to those of the VY2 plasma of 7 Cas. The ratios formed 
by these line strengths are more useful for these purposes 
than the fir triplets because the atomic level populations 
are found to be affected by UV photoexcitations by only 
«10% or less (Mauche 2002). 

In the competition between radiative and collisional de- 
excitations of the upper atomic levels, high densities mainly 
affect the lines arising from transitions from the lowest 
member of an excited manifold. For the resonance FeXVII 
lines we can observe these terms of these manifolds corre- 
spond to the 2p 5 3d and 2s 2 3s terms. The observable tran- 
sitions most affected are 15.26 A and 17.10 A, respectively. 
Both of these are too weak to be used in our analysis, so 
we must resort to alternative diagnostics. 

Fig. [7] exhibits the principal FeXVII ion lines visible in 
the 15-17 A region. The spectra shown are averages of fea- 
tures recorded in both RGS cameras and smoothed over two 
wavelength bins. The line ratio formed by 17.10 A (mag- 
netic quadrupolar transition) and 17.05 A is a density di- 
agnostic (Klapisch et al. 1978, Mauche et al. 2001). We 
estimate that in our spectra the ratio I(17.lA)/I(17.05A) 
lies in the range of 0.0-0.4 (with a value near zero most 
likely). This means that the electron density is of order 
10 14 cm -3 or higher, according to Mauche et al. (2001). 
For OBS2 the ratio I(16.8A)/I(17.0A) ps 1. From the work 
of Mauche et al. this ratio implies a density near 10 13-14 
cm -3 for both the high-flux and low-flux spectra. However, 
for OBS3 this same ratio has values ~1 for the high flux 
spectrum and decidedly less than 1 for the low flux spec- 
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trum. The last ratio indicates N e ;> 10 14 cm" 3 . (This same 
condition occurs for both the low and high spectra of OBS4, 
although it is not shown.) This difference suggests that the 
densities of the emitting plasma can sometimes fluctuate 1.0 
even on timescales of several minutes. Moreover, that the 
low-flux spectra constructed from any of the other obser- 
vations shows no comparable weakness of 16.78 A and that 0.5 
the higher density range is more typical. Altogether, the 
FeL ion lines suggest high densities for the &T2 plasma | 
component. These electron densities are centered near 10 14 § . 
cm" 3 and can fluctuate around a mean value of 10 13 " 14 M 
cm" 3 on either short or long timescales. The densities for ^ 
this plasma component are similar to those found for the _ 5 

high temperature flares from independent plasma cooling 
arguments (SRC). This high density constitutes one of the 
more surprising results we have found in this study. ^ 



7. Discussion 

7.1. Connections to the circumstellar environment 

Our analysis of CHARA MIRC and VEGA data resulted 
in similar disk solutions. There is also an indication that 
matter may be dispersed beyond the disk and is visible in 
the near-IR. As for the disk itself, its extent is larger in the 
H-band than the R-band (near Ha), which is consistent 
with the higher disk opacity expected in the H-band. 

In addition, only the MIRC ring model gave acceptable 
errors. Certainly, a ring model cannot be ruled out from our 
study. However, if future LPOI work demonstrates that the 
disk reaches the star, this result will have to be reconciled 
with the presence of spectral msf and brief UV continuum 
light dips These signal the presence of corotating struc- 
tures close to the Be star. It would seem that shear forces 
would destroy a high density structure in a transition zone 
between corotating regions and a Keplerian disk extend- 
ing inward that intrudes upon them. Thus, we expect the 
shears would force the inner disk transition to evolve to a 
ring structure. It may be important that Hony et al. (2000) 
found that the weak Brackett emission lines in the 7Cas 
spectrum are broadened to w 550/sini kms" 1 and gave ar- 
guments why the broadening is most likely kinetic. This 
velocity is significantly above the projected stellar rotation 
rate. This inference could be further evidence of corotating 
matter close to the star. In the event that the disk touches 
the star, the corotating clouds would have to lie elsewhere, 
e.g., outside the disk plane at intermediate latitudes. 

Fig. [5] gives a feel for the geometry of the system viewed 
from our vantage point. Models derived from LBOI are suf- 
ficient to describe the relative contributions of the 7Cas 
star/Be disk system. They indicate that its disk contributes 
22% to the V band brightness and only 3-4% in B, of 
the star-disk system (Stee & Bittar 2001). (Note that this 
prediction depends on disk conditions observed in 1993 
and on assumptions about the disk-photosphere transition.) 
The B — V color variations in FigQJi confirm these predic- 
tions. The Be disk is the only practical place in the system 
where secular and irregular continuum flux variations can 
originate in visible wavelengths. Moreover, the season-to- 
season averages shown in Fig.[TJi have a linear regression 
slope AB/AV = 0.66 ±0.02. This slope is consistent with 
color/brightness ratio changes observed when disks of clas- 
sical Be stars evolve in time (e.g., Hirata 1982). 



1.0 0.5 0.0 -0.5 -1.0 
Milliarcsecond 

Fig. 8. Artistic rendition depicting the IR H-band emission 
determined for 7Cas. We show the derived extent, elonga- 
tion, and position angle of the Gaussian disk. The inner 
region is optically thick. Neither partial obscuration of the 
star by the disk nor the star's gravity-darkening are in- 
cluded in the MIRC model. 



Table 5. Density columns (in 10 22 cm 2 ), Gaussian EW 
of Fe fluorescence, Abundance Anomalies (solar; Fe refers 
to K-shell) 



Date 


Binary 


n Hb 


EW(FeK) 


Fe 


Ne 


N 




Phase 




(mA) 








2001 


0.75 


10 


-19 


0.10 


~1 


~1 


2004 


0.27 


0.23 


-10 


0.12 


2.63 


3.96 


2010 


0.74-0.96 


36-74 


-35 to -50 


0.18 


1.80 


2.33 



"Fe" refers to Fe abundance from K-shell ion lines. 

Errors for [Fe], [N], and [Ne]: ±0.02, ±0.75, and ±0.28, resp. 



In t j6.2l we reported that during the 44 day interval of 
our 2010 observations the nn b column increased by about 
a factor of two. This rate of increase defines the column at 
which this range can change with time. Altogether a change 
of 300 x has been found between 2004 and 2010. As already 
stated, our observations covered the orbital interval be- 
tween inferior conjunction (Be star closest to the Sun) and 
RV maximum. The orbital phase (0.75) at which Chandra 
observed the large nu b in 200 f was almost identical to the 
phase of our first 20 f observation (0.74). Although it was 
still high in 2001, this column was smaller than in any 
of the 2010 observations. As the XMM-Newton sequence 
progressed toward quadrature, the njj b became higher yet. 
These values are summarized in Table 5. 

From this it is difficult to maintain that changes seen in 
riH b in 2010 can be attributed primarily to binary phase. 
If they do correlate with phase it must be in a complicated 
manner. Instead, as indicated in Table 5, nu b absorption 
correlates strongly with the reddening. We exhibit this re- 
lation in Figure^] This plot shows that the relation exists 
not just for the 2001 and 2004 observations but also for the 
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Fig. 9. The APT-measured B — V colors corresponding to 
the epochs of the six n# b absorption columns determined 
from soft X-ray flux attenuations. The n^ values increased 
monotonically during 2010. The errors in determining n# b 
are ±10%. 

four in 2010. Since we have associated reddening with disk 
development, this figure suggests that an association exists 
between X-ray absorption and disk conditions. 

We also find in Table 5 a suggestion of temporal vari- 
ance between the nu b column and the strengths of the Fe 
fluorescence feature, the K-shell-derived abundance of Fe, 
and of [N] and [Ne]. It appears that the status of the Be 
disk is related to several peculiar properties of the X-ray 
emitting plasma. 

Although it is by no means clear how these physical 
parameters are connected, one might begin to connect the 
dots by noting that the smaller than unity ratios of AB/AV 
observed by the APT are likely caused by the Be star 
injecting matter into its surroundings. In fact, variations 
in V and B of 0.015-0.025 magnitudes over several hours 
have been observed on a few occasions (SHV) . During these 
events the ratios of the AB/AV were typically 0.85-90, i.e., 
grayer than longer term variations as recorded in Fig.[TJi. 
Such events suggest a highly dynamic environment in some 
places over the Be star, perhaps as associated with failed 
local ejections or magnetic prominences. 

For more successful ejections to be associated with the 
riH b column, the ejections must: i) occur in lines of sight 
subtending a portion of the Be star, ii) be sustained over 
weeks or months, and iii) have a sufficient density that the 
emission measure (N e 2 integrated over volume) be visible 
in continuum light. For this study we cannot confirm the 
validity of these conditions from spatial models because 
our LBOI observations are made in optically thick wave- 
lengths and thus do not measure changes in volumetric den- 
sity. In consideration of ii) (sustainmcnt), observations of 
7 Cas in the IUE archives at MAST0 disclose that histori- 
cally this star often alternated between prolonged periods 
of strong and weak states of the wind emanating from the 
Be star, with short transitions of 1-3 months between them. 
Generally, though not always, the Ha emission EWs were 
found to correlate with the wind state for 7 Cas and other 
Be stars (e.g., Doazan et al. 1983, 1987). Therefore, it may 



5 "MAST" represents the Multi-Mission Archive at the Space 
Telescope Science Institute (STScI). Support for MAST is pro- 
vided by NASA Office of Space Science via grant NAS5-7584. 



be that the 2010 event we have reported for this star is 
one such "wind" episode. We also note that there is a dis- 
crepancy between the range found in the optical depths of 
the UV resonance lines and the larger range of viH b values. 
Lacking a new opportunity for simultaneous UV- X-ray ob- 
servations of 7 Cas, a resolution of this issue is not at hand. 

We make no assumption about what ejection mech- 
anism might be responsible for the mass loss or tossing 
episodes. As suggested in [21 a wind driven by radiative 
and centrifugal forces alone may play a role in this, but if 
so it is likely to be modified by additional processes that are 
not identified among classical Be stars. One attribute in the 
7 Cas wind not found in others is a "nickering" in the violet 
wing of the Si IV and C IV resonance lines on a timescale 
of 1\ hours (Cranmer, Smith, & Robinson 2000; "CSR"). 
This timescale is also seen in the X-ray flux of this star. 
CSR suggested that the rotation into view of X-ray active 
centers on the Be star increases the local wind ionization, 
thereby causing optical depth changes but not otherwise in- 
fluencing the wind structure or mass loss rate. Whatever the 
cause, disk diagnostics would be influenced by changes in 
the geometry and instantaneous rates of the outflows. Even 
so, these might not change the strength of the X-ray Fe-K 
fluorescence feature (formed by a wide range of Fe ions) 
and certainly not the Fe, N, and Ne abundance anomalies. 
Although none of this is understood, it is at least clear now 
that these abundances can vary on timescales of 1-2 years 
in the X-ray emitting plasma. This result cannot be under- 
stood as a result of recent nucleosynthetic processing in the 
stellar interiors. 

All three high resolution X-ray studies of this star have 
found it necessary to include two absorption columns into 
the X-ray fitting models. This raises the question of why it 
is that no matter what the riH h column density happens to 
be at any one time, it includes only the same fraction, 23- 
27%, of the X-ray emitting volumes. We can address this 
question by supposing that the X-rays are emitted on or 
very near the surface of the Be star from a confined range 
of latitudes. One can then imagine that as the flow escapes 
and crosses our line of sight it will clear all but the "lower" 
portion of the stellar disk. The lower region of the star 
covered by this outflow, perhaps 25%, would be the X-ray 
emission volume suffering the njj b absorption. In this pic- 
ture the fraction of the disk obscured by the column would 
be some function of the confinement of the outflow in lati- 
tude, the inclination of the star and Be disk, the oblateness 
of this very rapidly rotating star, and the separation of the 
inner edge of the Be disk from the star. With the right se- 
lection of parameters the areal coverage faction of the nn b 
column can remain approximately constant. 

Our picture must be confronted with additional optical 
photometry and X-ray spectra. The predictions for 7 Cas 
are clearly that sudden brightness increases and B — V red- 
denings should continue to be simultaneously accompanied 
by increased amounts of expelled particles. A portion of 
the released matter, perhaps this is an episodic wind, is di- 
rectly observable in the soft X-ray domain by its effect on 
this flux. If verified, this picture promises to be more evi- 
dence that some, and by implication all, of the hard X-ray 
flux of 7 Cas is created close to and in front of the Be star. 

7.2. A kTi-wind connection? 

Our analyses of the X-ray spectra have not been able to 
clearly identify the location of the formation of the sec- 
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ondary (fcTj-fcTa) plasma components. The FeL-shell ion 
line ratios yield high density for the fcT2 component, imply- 
ing that it is formed near the primary kT^ emission source. 
As to the cool fcTi component, L10 discussed the pro's and 
con's of whether it is associated with the standard wind 
of single hot stars. Their conclusion was the EMI of this 
component is larger than found for standard winds but its 
temperature is consistent with them. In this work we can 
take this conclusion one step further by pointing out that 
at least some of the slow l2-iofcev/Io.3-2fceV X-ray color 
variations shown in Fig. [3] arises from the soft energy band 
for which fcTi is an important contributor. Because rapid 
fluctuations have not been reported from X-rays originat- 
ing from shocks in B star winds, we believe this component 
does not arise primarily from winds: fcTi is somehow part 
and parcel of the X-ray generation processes in 7 Cas. 

8. Conclusions 

This study combines contemporaneous information ob- 
tained by a multidisclipinary campaign on a complicated 
7 Cas system (Be star + binary orbit + Be disk). We have 
refined its orbital ephemeris: it is in a nearly circular orbit 
with P or b ~ 203.53 days. We also confirmed the Meilland et 
al. (2007) result using new LBOI observations in Ha light 
that the Ha-emitting disk rotates at a rate consistent with 
Keplerian rotation. Its extent has been traced out to about 
one half the Be star's Roche lobe radius. We find from op- 
tical spectroscopy and imaging that there is no longer evi- 
dence of a visible one-armed disk structure. Moreover, the 
V an R emission components in the Ha line are weaker 
than before the year 2000. These changes are not atypical 
in the context of the behavior of Be disks. 

Yet in other respects 7 Cas's recent history is not typi- 
cal. The APT record shows an abnormal rate of brightening 
and reddening of 7 Cas, and there can be little doubt that 
this is caused by the interaction of newly ejected particles 
with the disk. Such events are often observed in monitoring 
Be stars, so this can have little to do with the star's binarity 
or the fact that it is emits copious X-rays. Fortuitously, the 
activity the APT recorded also coincided with the onset of 
our four XMM observations. As stated in M6.2I the attenua- 
tion of soft X-ray flux is caused by an increase in a column 
density n^ that covers 7 of X-ray emitting sites. The rate 
of increase in nn b during the 2010 observations is consistent 
with the changes in n^ measured in 2001 and 2004. These 
column densities in turn are correlated with the changes in 
the V magnitude and B — V color but are not related in 
a simple way to orbital phase. We argued that the logical 
place for the nn b to originate is from the time-dependent 
efflux from the Be star. We noted that these variations oc- 
curred without the additional variable column provided by 
a 1-armed disk structure. X-ray studies during the preces- 
sion of a future V/R cycle could test conclusively whether 
the iiH b variations are associated with the central Be star. 

Certain other aspects of the X-ray emitting plasma 
discussed in L10 are also clearer than before. The multi- 
thermal component nature of the sites seems relatively sta- 
ble. We stress that the soft X-ray flux from the hot fcT4 
component dominates the contributions of the cooler com- 
ponents. Even so, one can still see that the electron densi- 
ties of the warm plasma component are high, probably 10 13 
cm -3 or more. The cool X-ray plasma component clearly 
forms in different volumes but seems to have a higher emis- 



sion measure, even in 2004, than is expected for a typical 
wind of an early B star. In addition, we noted that the soft 
X-ray flux sometimes exhibits rapid aperiodic variations. 
This is not a known characteristic of Be star winds. 

From the correlation between color, and X-ray column 
density shown in Fig.O it is reasonable to conclude that 
j (all?) of the hot X-ray sites lie within a high density 
(>10 14 cm -3 ) medium like the photosphere of a star and 
the observer. In principle, an alternative interpretation is 
that the figure represents statistical flukes. A second al- 
ternative is that the column density represented by n# b is 
associated with accretion onto a dense column above a de- 
generate companion. There are difficulties associated with 
this picture, one of which is the faithful response of particle 
densities in the downstream column to an input variation 
far upstream in the Be wind. The picture would assume 
that the variations in the wind, the dense inner regions of 
the Be disk where the V flux is formed, and infall to a 
secondary some distance away, are all synchronized, with 
at most a short delay. Even assuming this chain of events 
offers no obstacles, such a picture would have to confront 
the correlations of rapid variations between X-ray and UV 
continuum and line fluxes expected to be visible only in the 
vicinity of the Be star. 

Aside from its correct interpretation, the reddening- 
column correlation promises for the first time to predict 
readily observable changes in X-ray properties of 7 Cas. It 
will be interesting to see with a future instrument if this 
relation extends to the column densities, abundances, and 
fluorescence strengths found for other 7 Cas members. 
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